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Impact phenomena of small clusters subject to thermal fluctuations are numerically 
investigated. From the molecular dynamics simulation for colliding two identical clusters, 
it is found that the restitution coefficient for head-on collisions can exceed unity when the 
colliding speed is smaller than the thermal velocity of one "atom" of the clusters. The 
averaged behavior can be understood by the quasi-static theory of impact phenomena. It is 
also confirmed that our result is governed by the fiuctuation theorem. 

§1. Introduction 

Inelastic collisions are the process that a part of the initial kinetic energy of 
colliding bodies is distributed into the internal degrees of freedom. Although it is 
believed that the relative speed of colliding two bodies decreases after a collision, the 
event that the relative rebound speed becomes larger than the relative colliding speed 
is not contradicted with the energy conservation law. Such an anomalous impact is 
prohibited by the second law of thermodynamics in which the energy does not come 
back to the macroscopic degrees of freedom once a part of energy is distributed into 
the microscopic ones. 

For classical colliding bodies at the initial temperature T, the second law may 
be written as 

^^V' + 0{T)>\^^V' (M) 

where /i is the reduced mass of the colliding bodies, V and V are respectively the 
relative colliding speed and the relative rebound speed. ^ The second term in the left 
hand side is negligible in the case of macroscopic bodies. Thus, in the case of head-on 
collisions of macroscopic bodies, the restitution coefficient, e = V /V, becomes less 
than unitj® while the restitution coefficient projected into the normal direction of 
the collision can exceed unity in the case of oblique collisions. 3''^' The second term 
in the left hand side of eq.f ll-ip . however, is not negligible for small systems. In this 
paper we investigate the statistical properties of "super-elastic" collisions for small 
colliding bodies where e is larger than unity in head-on collisions. 

"Super-elastic" collisions may be related to the fluctuation theorerrP''^''^®''^ 
in which the probability of the entropy production is related to the entropy absorp- 
tion as p^g^ JL^) = 6xp(rA), where P{S'^ = A) is the probability distribution of 
observing the time averaged entropy production S'^ during time interval r lies in the 
range A to A + dA. Thus, it might be important to study the relation between the 
fluctuation theorem and "super-elastic" collisions due to large thermal fluctuations. 

There are numerical and theoretical studies for coalescence, scattering, and frag- 
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mentation of colliding nanoclusters. HSJ'ttU'IIS Most of the low-speed collisions of 
nanoclusters cause coalescence. However, some stable clusters such as fullerene can 
keep their form in collisions. Therefore, we focus on the properties of small stable 
clusters in which the interaction between two clusters is dominated by the repulsive 
force. 

In this paper, we perform the molecular dynamics simulation of colliding repul- 
sive clusters to investigate the effect of thermal fluctuations. The organization of this 
paper is as follows. In the next section, we introduce our model. In fj3l we investigate 
the relation between the restitution coefficient and colliding speed and that between 
the compressive force and the deformation. We also compare our numerical results 
with the fluctuation theorem. Section [J] and [5] are devoted to the discussion and the 
conclusion of our results, respectively. 



§2. Model 



Let us introduce our numerical 
model. Our model is composed of two 
identical clusters. Each of them is 
spherically cut from a face-centered cu- 
bic (FCC) lattice and consisted of 682 
"atoms". The clusters have facets due 
to the small number of "atoms" (Fig. 
[1]). All the "atoms" in each cluster are 
bounded together by the Lennard-Jones 
potential U{rij) as 

(2-1) 

where rij is the distance between two 
"atoms", i and j, in one cluster, e is the energy constant and a is the lattice 
constant. When we regard the "atom" as argon, the values of the constants become 
e = 1.65 X lO^^^J and a = 3.4A, respectively.'^^ Henceforth, we label the upper and 
the lower clusters as cluster and cluster C', respectively. The interaction between 
the atom k on the lower surface of C" and the atom I on the upper surface of C' is 
assumed to be the repulsive potential R{rki) = Ae{a /vkiY'^ , where r^; is the distance 
between the atoms k and /. To reduce computational costs, we introduce the cut-off 
length dc of the Lennard-Jones interaction as = 2.5(t. 

The procedure of our simulation is as follows. As the initial condition of simula- 
tion, the centers of mass of C" and are placed along the z-axis with the separa- 
tion (Tc between them. The initial velocities of the "atoms" in both C" and obey 
Maxwell-Boltzmann distribution with the initial temperature T. The initial temper- 
ature is set to be T = O.Ole or T = 0.02e in most of our simulations. Sample average 
is taken over different sets of initial velocities governed by the Maxwell-Boltzmann 
velocity distribution for "atoms" . 



Z 




Fig. 1. Numerical model of colliding clus- 
ters. Each of them is composed of 682 
"atoms" which are bounded together by 
the Lennard-Jones potential. 
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To equilibrate the clusters, we adopt the velocity scaling methocP^''^ for 2000 
steps at the initial stage of simulations. We have checked the equilibration of the 
total energy in the initial relaxation process. After the equilibration, we give trans- 
lational velocities and the macroscopic rotations to C" and C' to make them collide 
against each other. The relative speed of impact ranges from V = 0.02y^e/m to 
V = 0.07-^e/m, which are less than the thermal velocity for one "atom" defined by 
y'r/m, where m is the mass of the "atom" . 

Numerical integration of the equation of motion for each atom is carried out by 
the second order symplectic integrator with the time step dt = 1.0 x 10~^cr/ y^e/m. 
The rate of energy conservation, \E{t) — Eq\/\Eq\, is kept within 10~^, where Eq is 
the initial energy of the system and E{t) is the energy at time t. 

We let the angle around 2— axis, 6^, he 9^ = when the two clusters are located 
mirror-symmetrically with respect to z = 0. In most of our simulation, we set 9^ at 
6^ = as the initial condition. Let us comment on the dependence of the relative 
angle 9^ on the numerical results. From our impact simulation for 9f = vri/lS (i = 
1, 9) at T = 0.02e we have confirmed that the initial orientation does not largely 
affect the restitution coefficient. 

§3. Results 
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Fig. 2. Relation between colliding speed and 
restitution coefficient. The solid and bro- 
ken lines are results of the quasi-static the- 
ory. 



h/o 0.36 




Fig. 



3. Relation between compressive force 
and deformation. Cross points and error 
bars are average and standard deviation of 
10 numerical results. Error bars are hardly 
seen due to small standard deviation. Solid 
line is the result of Hertzian contact theory 
for two spheres pressed each other. 



Figure [2] shows the relation between the relative speed of impact V/y^ e/m and 
the restitution coefficient e. Cross points and error bars are the average and the 
standard deviation of 100 samples for each value on x-axis. From this result, we 
confirm that the restitution coefficient e decreases with the increase of the colliding 
speed V/y^e/m. When the colliding speed is V = 0.02y'e/m at T = 0.02e, the 
average of e becomes 1.04 which is slightly larger than unity. It is interesting that 
our result can be fitted by the quasi-static theory of low-speed impacts 1 — e oc 
yi/^i6li,[T7},[i8} j£ ^j^g restitution coefficient at = is replaced by a constant larger 
than unity. Indeed, the solid and the broken lines in Fig. [1] are fitting curves 
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of e = «! — Q2 (^y/y^e/mj , where ai and 02 depend on material constants of 
colliding bodies. 

To remove the possibility of accidental agreement between the fitting curve and 
the data, we also check the validity of Hertzian contact theory in our system. Accord- 
ing to the theory, when the two elastic spheres are pressed each other, the relation 
between the deformation h and the compressive force P is described as h = DP^^^ 
and D = (3/2)2/3((l _ u"^) / Efl'^ , where u and E are Poisson's ratio and 
Young's modulus, respectively. US'® Poisson's ratio and Young's modulus of the 
model can be calculated as follows. 

Adopting these material constants, we confirm that our numerical result of com- 
pression of two clusters is well described by the Hertzian contact theory without any 
fitting parameters between P = 81.84e/o- and P = 136.4e/cr at T = 0.03e (Fig. [3]). 
For smaller compression finite deformation seems to remain. For larger compression 
we observe too large repulsion because of the existence of local plastic deformations. 
Thus, we conclude that the relation between the impact speed and the restitution 
coefficient is characterized by the quasi-static theory of impact processes. 

At last, we compare our numerical results with the fluctuation relation for impact 
phenomena which is a kind of the fluctuation theorem on the basis of the probability 
distribution of the macroscopic energy loss. In the case of impact phenomena, the 
fluctuation relation may be written as followsP 

exp(Vr(Xo,Xi)/r)P(VF(Xo,Xi)) = P(t^(Xi,Xo)), (3-1) 

where Xi{i = 0, 1) are the macroscopic variables at initial and final states, re- 
spectively, while Xi are the states obtained by reversing all the momenta in Xi. 
P{W{Xq, Xi)) is the probability distribution of W{Xq,Xi) which is the macroscopic 
energy loss during the transition from Xq to Xi defined as 



i=C^n'- 



\M{Vi-V'^)^\M-J^) 



+R{r) - R{r ). (3-2) 

Here, M is the total mass of one cluster. and Vi are respectively the angular 
velocities and the speed of the center of mass for the cluster i at initial state while 
and are respectively those at final state. / is moment of inertia of the cluster. 
R{r) represents the repulsion potential, where r is the distance between the centers 
of mass of the two clusters. 

In our simulation, we at first equilibrate the two clusters at T = 0.02e and collide 
them each other with the initial condition V = 0.02(e/m)^/^, = 10~^{e/m)^^'^ /a, 
and r = 12.02o". Here we do not use the initial macroscopic rotations induced by 
the thermal fluctuations with the average 9.5 x 10"''(e/m-)^/^/c7 and the standard 
deviation 4.9 x 10~''(e/m)^/^/cj as it is. Instead, for simplicity of our simulations, 
we give the initial uj whose value is approximately equal to \/Juj^ without taking 
into account the fluctuation of the initial rotations. We calculate W = W{Xq, Xi) 
from the initial and the final macroscopic energy by assuming the final state as the 
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state when the internal energy of the clusters keeps constant value in time. After 
the collision, we equilibrate the clusters at the initial temperature, and reverse the 
translational velocities to make them collide each other again. At the termination 
of the second collision, we calculate W = W{Xi, Xq). We obtain the probability 
distributions of W and W from 5000 samples. On the basis of the probability 
distributions, we investigate whether the relation ( 13-ip holds in this system. 




w W 

(a) (b) 

Fig. 4. (a) Numerical relation between In P{W)/P{W) and W for initial temperature T = 0.02e. 
Cross points are numerical results while the solid line is results from the fluctuation relation of 
inelastic impacts in ea.( I3-1|I . (b) The probability distribution of W and W obtained by our 
simulation. 

Figure SJa) shows the relation between W and In P{W)/P{W). The solid line 
corresponds to W/T at T = 0.02e. We find that our simulation data are nearly 
consistent with the line in the range of — O.Ole < W < 0.02e. Thus, the relation 
( 13-ip holds in the restricted range of W in our simulations. 

Here we should comment on the region in which our simulation results are con- 
sistent with the fluctuation relation. We show both P{W) and P{W) in Fig. [Hb). 
Both distributions are overlapped mostly in the range of — 0.02e < W < 0.02e. In 
the range of \W\ > 0.02e, the values of In P(W)/P{W) are calculated from few sam- 
ples. Thus, our data and the relation (|3-ip show good agreement only in the range 
—O.Ole < < O.Ole, while the discrepancy becomes large outside the region. Since 
the effect of rotations is small, W is approximately proportional to 1 — e^. Thus, 
events with the negative W almost correspond to "super-elastic" collisions. 

§4. Discussion 

Let us discuss our results. While our model mimics impact phenomena of small 
systems subject to large thermal fluctuations, we should address that our model may 
not be adequate for the description of most of realistic collisions of nanoclusters, 
where the adhesive interaction between clusters often prohibits the rebound in the 
low-speed impact. Thus, "atoms" in our model may be regarded as a coarse-grained 
units of crystalline solids and the dominant interaction among clusters is repulsive. 
We also expect that our model is an idealized collision of stable fulleren^ID' or charged 
clusters in which attractive interaction between clusters is weak.^^ As an additional 
remark, we should indicate that it may be diflicult to control the colliding speed 
and the initial rotation of the cluster in actual situations because the macroscopic 
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motion of one cluster is also affected by thermal fluctuations. 



§5. Conclusion 



In conclusion, we have performed molecular dynamics simulations to investigate 
the behaviors of colliding clusters and the relation between those and the fluctuation 
theorem. The results of our simulations have revealed that the relation between 
colliding speeds and the restitution coefficient can be described by the quasi-static 
theory for inelastic impacts. In addition, on the basis of the distribution function of 
macroscopic energy loss during collision, we have shown that our numerical results 
can be explained by the fluctuation relation for inelastic impacts. 
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